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Abstract 

In  this  paper  plane  elastic  curves  are  revisited  fronn  a  viewpoint 
that  emphasizes  curvature  properties  of  these  curves.  The  family  of 
elastic  curves  is  considered  in  dependence  of  a  tension  parameter  a  and 
the  squared  global  curvature  maximum  fcj^.  It  is  shown  that  for  any 
elastic  curve  is  bigger  than  the  tension  parameter  o.  A  curvature 
analysis  of  the  fundamental  forms  of  the  elastic  curves  is  presented. 
A  formula  is  established  that  gives  the  maximum  turning  angle  of  an 
elastica  as  a  function  depending  on  and  a.  Finally,  it  is  shown 
that  an  elastic  curve  can  be  represented  as  a  linear  combination  of 
its  curvature,  arc  length  and  energy  function  and  that  any  curve  with 
this  property  is  an  elastica. 


1  Introduction 


The  search  for  Smoothing  Algorithms  in  CAGD  has  induced  research  on 
curves  and  surfaces  which  minimize  functionals  with  geometrical  or  physical 
meaning.  The  classical  example  are  elastica  which  describe  the  shape  of 
elastic  materials.  Plane  elastica  can  be  defined  as  the  extremals  (critical 
points)  of  the  variational  problem 


+  O'  ds  — ♦ 


min 


where  k  denotes  the  curvature,  s  the  arc  length  of  a  plane  curve  x  and  a  is 
constant.  The  total  length  L  of  x  is  considered  to  be  variable.  The  set  of 
comparison  curves  of  the  problem  is  the  set  of  all  C®®  curves  in  the  plane 
with  fixed  endpoints  and  fixed  tangent  directions  at  these  points. 
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Lee  and  Forsythe  found  (see  [6])  that  the  curvature  function  of  an  elastica 
satisfies  the  differential  equation 

«"(5)  +  +  A'/c(s)  =  0,  KeR.  ( 1 ) 

We  begin  this  paper  with  a  short  derivation  of  this  result  and  a  proof  that  this 
equation  is  also  sufficient  for  a  curve  to  be  an  elastica.  We  will  then  show  that 
for  any  elastica  the  square  of  the  curvature  function  has  a  global  maximum 
which  is  bigger  than  the  tension  parameter  a.  The  fundamental  forms  of 
elastic  curves  first  described  by  Euler  are  characterized  by  the  values  of  the 
parameters  and  er.  A  curvature  analysis  of  these  curves  is  presented. 

In  section  3  various  formulas  expressing  properties  of  elastica  in  terms  of 
and  a  are  given.  This  includes  a  formula  that  gives  the  maximum  turning 
angle  of  an  elastica  as  a  function  depending  on  and  a.  The  explicit 
formula  for  the  maximum  turning  angle  allows  to  determine  an  upper  bound 
for  k'^  in  the  case  that  an  inflectional  elastica  has  a  turning  angle  bigger 
than  It.  As  the  main  result  of  this  section  we  show  that  a  plane  elastica  has  a 
representation  as  a  linear  combination  of  its  curvature,  arc  length  and  energy 
function  E{s)  =  Jq  «^(s)  ds.  This  property  characterizes  the  class  of  plane 
elastica  and  can  be  used  for  speeding  up  the  computation  of  interpolating 
elastica. 


2  The  curvature  of  plane  elastica 

Given  two  points  P,Q  €  R^  and  two  unit  vectors  V  €  TpR^,W  6  7qR* 
then  M  denotes  the  set 

M  :=  {x  :  [0,  L]  :  I  €  R+,x  e  C~[0,  L],  |x'(s)|  =  1  for  s  €  [0,1], 

x(0)  =  P,  x{L)  =  (?,  x'(0)  =  V,  x'(I)  =  W}. 

We  consider  the  problem  of  minimizing  the  functional 

£  +  <tL  :=  f  K^(s) -I- <7  ds 

among  all  curves  of  M  where  k  denotes  the  curvature  and  s  the  arc  length 
parameter  of  a  plane  curve  x.  The  integral  E  =  Jq  ds  is  called  the  en¬ 
ergy  of  X  on  [0,  L].  a  denotes  a  constant  which  we  call  the  tension  parameter. 
The  total  length  L  is  considered  to  be  variable.  As  the  functional  considered 
is  translation  and  rotation  invariant  we  assume  in  this  section  that  P  =  (0,0) 
and  V  =  (1,0). 
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For  X  €  M  the  tangent  vector  F  of  x  is  given  by 
r(a)  =  (cos(«I'(a)),  sin{*{5))) 

where  the  function  with  :=  k{s)  ds  gives  the  turning  angle  of  x. 
Using  this  notation  the  variational  problem  can  be  written  as 

fL 

min  /  +  (t  ds 

*€C<»10.L]Jo 

under  the  constraint 

I 

where  the  admissible  functions  €  C®®fO,  L\  are  subject  to  the  constraints: 

^(0)  =  0,  =  with  (cos0,  sin0)  =  W. 

According  to  Lagrange’s  multiplier  rule  the  differential  equation  of  this 
problem  is  the  Euler  equation  for  the  integrand 

Fi'b,  4^',  a)  =  +  O’  +  A  cos  ^>(s)  -f  (xsm  ^(s), 


for  some  constants  A,/i  €  R,  i.e. 

=  — ~sin'P(s)  +  ^cos^(s). 

(see  e.g.  in  [3], [6]).  By  introducing  the  constants 

A  =  2a  cos  4>  /i  =  2a  sin 
we  rewrite  the  Euler  equation  in  the  form 

=  k'  =  — asin(\I'  —  d>). 

Multiplying  (3)  by  2’!''  and  integrating  yields 

=  2a  cos(^'  —  ^)  +  A 


(2) 


(3) 


(4) 


TTAlV'iV'  ^  aaa  w«*aa  w  »  J  a-icw 

to  define  elastic  curves  in  the  plane  in  the  classic  literature  (see  {7]). 

In  order  to  determine  the  constant  A  in  terms  of  the  tension  par 
a,  we  consider  the  boundary  condition 

F(»(i),  rW,  l)  -  4;f (»(*),  ♦'(»),  =  0 


Per 

•I 
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Special 


that  must  be  satisfied  by  the  extremal.  This  condition  is  implied  by  the  fact 
that  the  total  length  L  of  the  curve  is  variable  in  the  variation  (see  e.g.  [2], 
p.  571).  Thus, 

=.  k’(I)  =  A  cos  ^{L)  +  /I  sin  4'(L)  +  a.  (5) 

Comparing  (5)  with  (4)  shows  that 

A  =  a. 

Therefore  we  give  the  following  definition  of  an  elastic  curve  under  tension. 

Definition  1  An  arc  length  parametrized  plane  curve  x  with  curvature  func¬ 
tion  K  is  called  elastica  (or  elastic  curve)  with  tension  parameter  a,  if  for 
some  a,  6  R 

=  2acos(4'  —  ^)  +  <T  (6) 

where  denotes  the  function  4'(s)  :=  K(s)ds. 

Lee  and  Forsythe  showed  in  [6]  that  the  Euler  equation  (2)  implies  the 
differential  equation  (1)  for  the  curvature  function  k.  We  will  now  show  that 
(1)  is  in  fact  equivalent  to  (6). 

Theorem  2  If  k  £  C^(R)  and  ’I'(s)  :=  /q  K{s)ds,  then 

—  2a  cos('f  —  d>)  + 

holds  for  some  constants  a,  9  and  o  if  and  only  if 

k"  +  -  ioK  =  0.  (7) 

Proof:  (i)  Differentiating  (6)  yields  the  Euler  equation  in  the  form  of 
(3).  Differentiating  the  Euler  equation  gives 

k"  =  —a  cos(’J'  —  (I>)k. 

Substituting  the  term  cos(^'  —  ip)  in  the  above  equation  according  to  (6)  we 
obtain  (7). 

(ii)  Differential  equation  (7)  can  be  integrated  to  the  first  oder  differential 
equation 

(«r  =  C-(l/4)<r’-(l/4)(«'-<T)’.  (8) 

Note,  that  for  a  real  solution  of  (8)  it  is  necessary  that 

a^:=C-(l/4)<7^>0 


4 


and 


4a2  - 


The  function  B  defined  by 

B  :=  (--l)’*arccos( 


ic^-a 

2|al 


) 


with  n  =  1  for  k'  >  0  and  n  =  2  for  k'  <  0  obviously  obeys  the  relation 

ic*  =  2(atcos5  +  <T.  (9) 

It  remains  to  show  that  B  is  an  integral  of  k. 

(8)  together  with  the  definition  of  and  (9)  implies 

(/c')^  =  a^sin^0. 

Since  sin  B  is  positive  (negative)  if  k*  is  negative  (positive),  we  obtain 

k'=— |a|sin^.  (10) 

Differentiating  (9)  and  substituting  k'  according  to  (10)  yields 

&  =  K. 


Since  'P  and  B  are  both  integrals  of  k  there  is  a  constant  (j>  such  that 
'i!  =  B-\-<i>.  For  ’i' 

=  2|al  cos(4'  ~  4>)  +  a 

holds  because  of  (9). 

□ 


The  squared  curvature  function  k’  of  an  elastic  curve  has  a  global  maxi¬ 
mum,  even  if  k*  is  extended  to  the  whole  real  line. 

Lemma  3  If  k  is  a  solution  of  (7)  on  R,  then  has  a  global  maximum. 

Proof:  Since  the  function  (l/2)(#c®  —  <tk)  is  of  the  class  C*(R)  the  solu¬ 
tions  of  (7)  can  be  extended  to  the  whole  real  line. 

To  show  that  has  a  global  maximum  we  observe  first  that  (3)  implies 
the  existence  of  a  local  extremum  of  «.  This  is  because  the  assumption 
k'(s)  =  '^"(s)  ^  0  for  all  s  €  R  means  that  is  convex  or  concave  and 
therefore  unbounded  while  according  to  (6)  k*  has  no  zeros  only  if  ^  is 
bounded.  Furthermore  from  (3)  and  (4)  it  is  obvious  that  any  local  extremum 
of  K  is  a  global  extremum  of  k?. 
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We  assume  now  that  has  no  global  maximum  on  R.  In  this  situation 
K*  has  exactly  one  zero  on  R  and  takes  its  global  minimum  in  this 
point.  Note  that  k  is  non-zero  for  any  s  ^  5^^  because  a  zero  of  k  at  a 
point  s  ^  s„,>,  would  imply  that  k(s„,„)  =  0  and  K'(i)  —  0  for  some  point  s 
between  s  and  Smin-  Therefore  (/c^)'  has  no  zero  besides  s„i„  and  k}  is  strictly 
monotone  increasing  on  the  right  of  Smin-  Hence  k  is  monotone  increasing 
resp.  decreasing  on  the  right  of  Smm  if  «  has  positive  resp.  negative  values 
on  the  right  of  Sm.n-  The  formula 

'I>(s)=  /  #c{s)ds-|- 

yields  that  ^  is  in  any  case  unbounded.  This  is  a  contradiction  to  the  as¬ 
sumption  because  (4)  implies  the  existence  of  global  maxima  for  if  is 
unbounded. 

□ 


We  now  give  the  relation  between  the  global  maximum  of  k?  and  the 
tension  parameter  a  and  express  k  in  terms  of  elliptic  functions. 

Theorem  4  Lei  k  6  C'^(R)  be  a  solution  of  the  differential  equation  (7)  with 
a  global  maximum  0  of  on  R.  Then  the  following  statements  hold: 
(0 

(ii)  K  has  a  zero  if  and  only  if  >  2<7.  In  this  case  k  is  given  by 


k{s)  =  K„  cn(^«  ~<r)/2(s  -  s,„)  |  k^) 


(11) 


with  the  parameter 


k^  = 


2(r  —  kI^<  K?  <  Kl^  for  <  2<7.  In  this  case  k  is  given  by 
k(s)  =  K^dniKnis  -  Sm)/2  I  p). 


(12) 


Proof;  As  a  global  maximum  point  s^  of  k*  is  also  a  zero  of  k'  (8)  implies 


Hence  (8)  takes  the  form 


/-I  _  1  1  —.,2 

4  m  2  *" 


1 


(«')’  =  4('«m  -  +  «m  -  2<t). 


(13) 


m 


4 
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As  all  quantities  in  (13)  are  real  (i)  follows  from  the  fact  that  the  term 
(/c^  —  K^)  is  always  non-negative  and  therefore  the  term  (k*  +  ~  2a)  has 

also  to  be  non-negative  for  any  k^. 

Using  the  same  argument  we  observe  that  k}{s)  =  0  for  some  s  implies 
that  Kl^  >  2<t  while  for  <  2a  #c^  has  to  be  greater  or  equal  2a  — 

To  express  k  in  terms  of  elliptic  functions  we  precede  as  follows.  In  the 
case  that  >  2a  we  substitute  —  k})Ik\  in  (13)  and  obtain  for  z 

the  differential  equation 

=  (14) 

where 


(14)  is  the  differential  equation  of  Jacobi’s  function  sn=sn(u)  for  the  argu¬ 
ment  u  =  —  a)s  (see  [l],p.  114).  We  therefore  obtain 

2^(s)  =  sn^(^(/c^  -<y)/2(5  -  5m)  I  ^*)- 

The  relation  =  /c^(l  —  z^)  then  gives 

=  (15) 

(see  [8],p.l6).  Since  k  is  differentiable,  (15)  implies  (11). 

In  the  case  that  <  2a  we  substitute  z^  =  («^  -  k^)/2(#cJ^  -  a)  in  (14) 
and  obtain  for  z  the  differential  equation 


where 


rt  _  1  _  2(4  -  <^) 

•  ««  • 


Hence 


The  relation 


r*(s)  =  sn^(Km(s  -  s„)/2  |  /^). 

=  «m(l  - 


yields 

=  «mdn^(«m(5  -  Sm)/2  ]  /^) 

(see  [8],p.l6)  which  implies  (13)  again  because  of  the  differentiability  of  k. 
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Theorem  2  implies  that  the  curvature  function  of  an  elastic  curve  extends 
to  a  periodic  function  on  R.  In  the  case  that  >  2<t  this  periodic  function  is 
symmetric  with  respect  to  any  zero  of  its  derivative  and  antisymmetric  with 
respect  to  any  zero.  According  to  Love  [7]  this  situation  is  called  inflectional 
because  the  extension  of  x  has  turning  points.  In  the  case  that  <  2a  the 
extension  of  k  has  no  zeros  but  is  still  symmetric  with  respect  to  any  zero 
of  its  derivative.  The  local  extrema  in  this  case  are  ic„  and  ±(2a  —  k^). 
The  situation  is  illustrated  in  figure  1  •  figure  8  where  for  «„  =  1  various 
curvature  functions  and  the  corresponding  elastica  are  shown.  If  only  non¬ 
negative  tension  values  are  considered  a  is  an  element  of  [0, 1].  With  a 
increasing  from  0  to  1  the  curvature  function  changes  continuously  from  a 
lemniscate  function  (a  =  0)  to  a  constant  (a  =  1).  A  classification  of  the 
different  forms  of  elastica  has  beei  given  by  Euler  (see  [5],  [7]). 
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3  A  characterizing  representation  theorem 

First  we  summarize  the  main  formulas  for  elastic  curves.  In  order  to  avoid 
repetitions  we  use  the  following  convention:  if  a  curve  x  :  [0,  L]  — ♦  R’  is  an 
elaistica  ,  then  k  denotes  the  (analytic)  curvature  function  of  i,  k  denotes  the 
extension  of  k  on  R,  the  global  maximum  of  and  s„  is  a  number  with 
Furthermore  ♦  :  fO,  L]  — »  R  is  the  function  4'(s)  :=  K{s)ds, 
’i'  denotes  the  extension  of  on  R  and  (p  :=  The  energy  E  of  i  is 

the  function  E{s)  =  jQK^(s)ds.  For  convenience  we  also  use  the  notation 
=  (cost/?, sin t/>).  t?  is  the  angle  in  [0,2x{  with  x'(0)  =  C(i?). 

Theorem  5  For  a  plane  elastic  curve  x  with  tension  parameter  a  the  fol¬ 
lowing  relations  hold: 

(i)  k'{s)  =  -5(k^  -  <r)  sin('I'  -  ip), 

(ii)  K^{s)  =  (k^  -  <t)cos('I'(s)  -  s?)  + 

(Hi)  k(s)  -  k(0)  =  -5(«^  -ff)  <  C'(ip  +  t>),x(5)  -  x(0)  >. 

(iv)  E{s)  =  {kI,  -a)<  C{ip  +  t?),x(s)  -  x(0)  >  +crs, 

Proof:  Without  loss  of  generality  we  can  assume  that  a  in  (4)  is  positive. 
From  lemma  3  we  know  that  k}  has  a  global  maximum  Since  for  positive 
a  the  maximum  occurs  iff  cos('I'(s,„)  —  <6)  =  1,  it  follows  that 

2a  =  Kl,-cr 


(j)  =  4'(s,„)  +  2kn. 

Therefore  (3)  and  (4)  imply  (i)  and  (ii).  To  obtain  (iii)  we  rewrite  (i)  as 

k'{s)  =  ~  <^)(sin(^'  +  t?)cos(v>  +  i?)  —  cos(4'  -f  sin(v?  +  t?)) 

and  integrate  using 

x'(s)  =  (cos(4'(s)),sin(fl'(s))). 

To  verify  (iv)  we  write  (ii)  in  the  form 

Ac^(s)  =  (k^  —  cr)(cos('l'  +  i?)cos(<,?  +  t?)  +  sin(^'  +  i?)sin((^  +  i?))  +  (T 
and  integrate. 

□ 
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Since  the  elliptic  function  cn  in  Theorem  3  has  the  symmetry  properties 
of  a  sine  wave  while  dn  is  positive,  the  turning  angle  of  an  elastic  curve  is 
bounded  in  the  inflectional  case  a  <  and  unbounded  otherwise.  Formula 
(ii)  of  Theorem  3  can  be  used  to  determine  the  maximum  turning  angle  of 
an  inflectional  elastic  curve. 


Corollary  6  The  maximum  turning  angle  •=  of  an 

elastica  x  with  period  T  and  tension  parameter  ff  <  is  given  by 

'I' mar  =  2  arccos{ — (16) 


If  a  >  ^  0  then 

for  any  t  6  R. 


l^(t  +  T)  -  =  2ir 


Proof:  Let  sq  be  a  zero  of  k.  Then  it  follows  from  the  symmetry  prop¬ 
erties  of  k  in  the  inflectional  case  that 

'I'mar  =  2|  *1  =  2|a'(So  +  r/4)  -  <'(o)l. 

J»0 


With  (ii)  we  get 

'I'max  =  2|  ATCcosi^ — -)  +  -  arccos(—  -^)  -  «^| 

and  therefore  (16). 

In  the  non-inflectional  case  k  is  a  positive  periodic  with  the  period  T  = 
\Kj\Km\  where  K  denotes  the  complete  elliptic  integral  of  the  first  kind.  For 
i-=  Sm  one  obtains 

|’^(s„,  +  T)  -  4'(s„)|  =  4  /  dn(ul/^)<iu  =  4arcsin(sn(/C))  =  2ir. 

Jo 

□ 


For  an  inflectional  elastica  i  formula  (16)  implies  an  upper  bound  for 
if  the  absolute  value  of  the  oriented  angle  between  i'(0)  and  r'(L)  is  bigger 
than  It. 


Corollary  7  For  an  inflectional  elastica  x  with  tension  parameter  cr  and 
:=  \'9{L)\  >  7t  : 


2  < 


<  1  - 


1 


cos(V»/2) 


4 
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•> 


) 


Proof:  It  follows  from  (16)  that  a  turning  angle  bigger  than  n  can  only 
happen  for  positive  <7.  The  inflectional  nature  of  x  implies  therefore  the  left 
inequality.  The  right  inequality  follows  from  the  fact  that  tp  has  to  be  less 
or  equal  ’i'mo*  which  is  given  by  (16). 

O 


As  the  main  result  of  this  paragraph  we  show  that  an  elaistica  can  be 
represented  as  a  linear  combination  of  its  curvature,  energy  and  arc  length. 

Theorem  8  If  x  is  an  elastica  with  tension  parameter  a  ^  then 

U  cos(9  +  tf )  sm(v,  +  fl)jl  E(s)-„s 
where  cos^p  —  (kq  —  <t)/(k^  —  a)  and  sinv?  =  2kqI{k^  —  <r). 

Proof:  A  plane  curve  x  €  C^’fO,  L\  has  a  representation 


X(s)  —  X(0)  =  /  C(4^(s)  —  (,?  +  (,9  +  I?)  ds 

Jo 

=  C{ip-\-^)j  cos(4'(s)  —  i,9)ds 
+  C sin{^ {s)  -  tp)ds  . 

Applying  (iii)  and  (iv)  of  Theorem  3  yields 


x(s)  -  x(0)  =  — 1— (C(.,9  +  0)iE{s)  -  cTs)  -  +  t?)2(«(s)  -  /c(0))) 

which  is  equivalent  to  (17). 

The  formulas  for  ip  follow  from  (i)  and  (ii)  for  s  =  0. 


□ 


Formula  (17)  provides  an  explicit  representation  of  an  elastica  in  terms 
of  its  curvature  function  if  ^  a.  Note  that  in  the  case  =  <t  the  elastic 
curve  is  a  circle  of  radius  l/|Km|  as  shown  in  figure  8. 

Finally  we  show  that  elastica  are  essentially  the  only  curves  in  the  plane 
that  have  a  representation  of  the  form  (17). 

Theorem  9  Let  k  be  an  arbitrary  C^(R)  function  with  a  global  maximum 
of  K^,  a  a  real  number  smaller  than  and  E{s)  :=  Jq  K^{s)ds.  An  arc 
length  parametrized  curve  x  where  z(s)  is  given  by  (17)  is  on  elastica  with 
curvature  function  a  and  tension  parameter  <r. 
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Proof:  We  assume  the  case  ^  =  »/2  which  can  alw'ays  be  achieved 

by  applying  a  rotation  to  x.  Then 

x'  =  -i— (2k',  k’  -  o),  x"  =  -i— (2k",  2kk'). 

m  ffl 

X  is  arc  length  parametrized  if  and  only  if  x'  =  1,  i.e. 

i2Kr  =  iKl-ar-iK^-a)\  (18) 

Differentiating  (18)  yields 


2k' k"  +  k'k(k^  —  <t)  =  0. 

Therefore  either  k  satisfies  (7)  or  k  is  constant.  For  an  arc  length 
parametrized  curve  the  determinant  [x',x"]  is  curvature.  Here  we  have 

-  <^r 

If  K  is  constant  the  curvature  of  x  is  zero,  hence  x  is  a  trivial  elastic  curve. 
If  K  satisfies  (7)  we  substitute  in  (19)  k"  according  to  (7)  and  k'*  according 
to  (18).  This  yields 

[X',X'1  =  K. 

□ 

The  representation  formula  (17)  of  an  elastic  curve  is  extremely  useful  for 
the  computation  of  interpolating  elastica.  One  reason  for  this  is  that  (17) 
involves  no  trigonometric  functions  and  fewer  integrations  than  the  standard 
representation  based  on  the  formula  x'  =  (cos ’i',sin 'I').  Therefore  (17)  is 
less  expensive  to  evaluate.  Furthermore  (17)  can  be  used  to  find  piecewise 
polynomial  approximations  of  elastic  curves  based  on  a  spline  approximation 
of  the  curvature  function.  The  author  has  established  polynomial  splines 
which  approximate  the  curvature  functions  of  plane  elastica  with  high  ac¬ 
curacy.  Using  these  piecewise  polynomial  curvature  functions  together  with 
the  new  representation  (17)  one  obtains  polynomial  spline  approximations  of 
the  elastica  itself.  These  approximations  will  be  discussed  in  [4]. 
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